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Abstract 

We review some recent results concerning the evolution of a vortex filament and 
its relation to the cubic non-linear Schrodinger equation. Selfsimilar solutions and 
questions related to their stability are studied. 
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1 Introduction 



1.1 The binormal flow 

This paper is mainly concerned with the so called Localized Induction approximation 
for the evolution of a vortex filament. In geometrical terms it is given by the binormal 
flow (BF): consider the 3d curve x( s ,t) £ ^ 3 with t £ K denoting time and s£l 
the arclength parameter, then the flow is given by 

Xt Xs A Xss /-i \ 

x (s,0)= Xo (s). [1) 

Denoting c the curvature and b the binormal vector we get from the Frenet equations 

' Xt = cb 

X (s,0) = Xo(s). [) 

This flow appeared for the first time in a work in 1906 by L.S Da Rios [DaRj as a 
model which describes the flow of a vortex filament. That is to a say a vortex tube of 
infinitesimal cross section which propagates according to Euler equations. Let us call 
x(s, t) a curve in 1R 3 which can be either closed, and then we will consider s £ [0, 2n] 
and 27r-periodic functions, or which extends to infinity in both directions. In that case 
s £ R. Let us consider for simplicity of the exposition this second possibility. 

Assume that the vorticity w of a given fluid of velocity u is supported along the 
filament described by the curve x(s,t) . By that we mean w is the singular vectorial 
measure 

w = TTds, 

where T stands for the circulation, 

T = Xs, 

s is arclength parametrization and T is constant along the filament. The velocity can 
be recovered from the vorticity through the Biot-Savart integral 

u(P) = f f°° { X ( S) ~L A T(s)ds. (3) 
4vr \x(s) - P\ d 

Typical examples of vortex structures are: 

1. The straight line vortex, x( s ) = (0,0, s). The velocity at a point P = (x,y,z) 
with z 7^ can be easily computed from ([3]) 

ds „ N C 



u(x, y, z) = -J-H/, x, 0) ! — — — — T = {-y, x, 0) 

47r J-oo x 2 + v 2 + z - s) 2 )2 



2 i „2 ' 



(x 2 + y 2 + (z - s) 2 )2 ** + V 

Hence the velocity around the filament is just a rotation of speed a constant times 
the inverse of the distance to the filament. As a consequence the filament remains 
stationary. 

2. The vortex ring, x( s >0) = (cos s, sin s, 0). 
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3. The helical vortex, x( s >0) = (cos s/v2, sins/ v2, /\/2),. -see |LFj and [Har| 

The computation of the velocity field outside of the filament is much more delicate 
in the last two examples and in fact, it is better to compute the so called vector 
potential using cylindric variables. We refer the reader to [AKQ], sections 2.5 
and 2.6 for the details. 

Notice that the Biot-Savart integral is perfectly defined for points P outside of the 
filament. The idea is to compute u(P) as P approaches the filament to eventually 
identify that limit with the velocity of the filament. Also observe that the term inside 
of the integral in (|3j) has enough decay at infinity. So in order to understand the 
leading behaviour we can just look at a neighbourhood of x( s o)> the closest point in 
the filament to P. A quick look at the integral is enough to show that the limit 
diverges as the inverse of the distance. In fact we see that the first approximation is 

lx(«o) - p\ 

This term behaves as in the case of example 1 of the straight filament. As we said 
the overall effect on the movement of the filament is null. Doing the Taylor expansion 
around that point, calling e the distance of P to the filament and removing the first 
term we get (see for example [Saj ) 

u(P) ~ —cT log e Xs(so) A Xss{sq) H 

Hence we see we have a logarithmic dependence on the distance to the filament on 
the direction of 

Xs f\Xss = b(s ). 
After a renormalization of time we get 

Xt = Xsf\Xss- (4) 

This computation is extremely crude and at best indicates that close to the fila- 
ment the leading term is in the direction of the binormal. The general consense (see 
[Sa| ) about this approximation also known as local induction approximation, is that 
it is specially suitable for nummerical integration but it has serious flaws and is of 
limited physical interest. In any case it would be extremely useful to find a way of 
desingularizing the integral in ([3]). 

Our next step is to establish the connection between the binormal flow and non- 
linear Schrodinger equation. A first hint of such a connection comes by identifying the 
normal plane to the curve x( s ) a t a given point with the complex plane. Then if n 
and b are the normal and the binormal vector respectively we can write the binormal 
as 

b = in. 

Hence the binormal flow can be understood as 

1 

- xt = cn. 

i 

Finally if we remember that the "curvature flow" has as a linear approximation the 
heat flow, we can "guess" that the underlying linear equation for will be Schrodinger 
equation. In fact something even stronger is true. 
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H. Hasimoto in |Hasj observed that if we define the complex valued function u as 

u(s,t) = c{s,t)e i ^^ s '^ ds \ 

where r stands for the torsion of the curve and c for the the curvature, and if x{ s ^) 
solves ([!]) then u solves 

iu t + u ss + ^(\u\ 2 + A(t))u = 0, (5) 

with 



A(t)-[2 Css + c 2 )(0,t). (G) 



Conversely, let A(t) be a time depending function, and u a regular non vanishing 
solution of ([5]). If we write 

u(s,t) = c{s,t)e i ^ s > t \ 
by defining r(s,t) = cp s (s,t), we have 

u(s,t) = ctMK^ r(s ''' )ds ' +<K °'* ) . 
We obtain then a filament function 

u(s,t) = c{s,t)e i ^< s '> t W } 

satisfying ([5]) with A{t) replaced by A(t) + (/>t(0, t). The real functions c and r verify 
then the intrinsic equations and the identity © with A(t) replaced by A{t) + 4>t(0,t). 

Notice that taking u = uexp ^— | J* A(t')dt'^j we get 

1. l2 

m t + tt ss + -\u\ u = 0, (7) 

which is the well known cubic non linear Schrodinger equation. 

In doing this H. Hasimoto uses the Frenet frame (T, n, b) and the corresponding 
Frenet equations 

T s = cn, 

n s = -cT +rb, (8) 

b s = — rn . 

There is no reason why to do this and more general frames (T, e\ , e^) can be used so 
that the assumption c > is not necessary, see for example |Ko| . The final conclusion 
is still the same and equation ([5]) is also obtained. 

Also it is important to observe that if in (J4|) we differentiate with respect to s in 
both sides we get, since Xs = T, 

T t = Xs /\Xsss = TAT ss . (9) 

Hence 

d . .o d 

-fa\Xs\ = 2^ Xs ■ Xs = 2(x s A Xss) s ■ Xs = 0, 

and the arclength parametrization is preserved. This is generally admitted as one of 
the flaws of BF as an approximation of ([3j) . 
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Notice that |T| 2 = 1 and T A T ss = JD S T S where J is the complex structure on 
the sphere and D s the covariant derivative. Therefore we can write ([5|) as 

T t = JD S T S 

which is nothing but the Schrodinger map onto the sphere in the one dimensional case. 
For higher dimensions we get 

T t = T A AT (10) 

which also appeared in ferromagnetism where is referred as the Lipschitz-Landau equa- 
tion. 

Let us finish this short introduction to the binormal flow observing that Hasimoto 
transformation allows us to exhibit some particular examples of (T5]): 

1. The straight line is just c(s, t) = 0, A(t) = 0. 

2. The circle is just c(s, t) = 1 A(t) = —1. 

3. The helix is just 

u(s,t) = e- itN2 e iNs A(t) = -l. 

4. The travelling wave solutions 

u( s , t) = e- itN2 e iNs Q(s - 2Nt), A(t) = -1, 

with Q(s) = sech(s). This last example is particularly interesting because 
suggests the possibility of the existence of travelling wave solutions along vortex 
filaments. In fact a few years after H. Hasimoto's paper, E.J. Hopfinger and F.K. 
Browand [HB] (see also [MHR ) gave experimental evidence of such a possibility. 

It is also interesting to recall some results by R. Klein, A. Majda and K. Damodaran 
[KMD] about the interaction of more than one filament. They propose a mathematical 
model based on a system of coupling NLS. This model is obtained as an equilibrium 
between the interaction of a filament with itself, that is described just by the linear 
part of BF, and with the others, that is measured by a sum of culombic type potentials, 
sec [MB] for the details. This system allows for helixes as solutions. In the case of two 
or three helixes with the same vorticity and at the same distance the configuration is 
stable, see [KPV1| . There is some experimental evidence that the case of four helixes 
could be unstable, see [AKOSj . 

Finally recall that the 1? norm is conserved in NLS. In our setting this becomes 

/oo /-oo 
c 2 (s,t)ds = / \u(s,t)\ 2 ds, (11) 
-oo J —oo 

which can be understood as the kinetic energy. Notice that in the case of the helix we 
obtain a solution of infinite energy eventhough it has finite L 2 -norm when consider as 
a periodic solution. As a filament one has to look at the complete helix because any 
vortex filament either ends in the boundary which confines the fluid or it has to go to 
infinity. This suggests as a natural question to consider solutions of (llip with infinite 
energy. 
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1.2 The mKdV flow 



The connection between NLS, which is a completely integrable system, and a geometric 
flow raised the natural question about the existence of more examples with similar 
connections. G. Lamb in the 70 's ([La]) was the first one in extending Hasimoto's 
idea, obtaining a remarkable generalization which include NLS and mKdV (even as 
complex valued functions) as possible limits. It is what is known as Hirota equation. Y. 
Fukumoto and T. Miyazaki [FM| obtained the same flow when in da Rios approximation 
is included the existence of an axial flow along the filament. So it is still a 3d-flow. 

Almost at the same time R.E. Goldstein and D.M. Petrich [GPlj proposed the 
mKdV hierarchy as the PDE of the curvature k(s, t) of plane curve flows which preserve 
the enclosed area. In a later paper and in a complete formal fashion they relate it to 
the dynamics of a vortex patch in the plane. Recall that for plane fluids the vorticity 
is a scalar, see [GP2]. A vortex patch describes the situation of a compactly supported 
vorticity with a jump discontinuity through the boundary. Assume for simplicity that 
the vorticity is a constant T in the patch f2 which is interior to the boundary d£l . 
Then N.J. Zabusky, M.H. Hughes and K.V. Roberts [ZHR] obtain an equation for the 
boundary which is given by a singular integral coming from Biot-Savart law. After 
doing a formal Taylor expansion in this integral (see [GP2]) the first term is just a 
translation along the boundary. The second term is 

z t = — z.sss + z ss z t = k s n — —k^T. 

If we call k the curvature one gets 

3 2 

kt = k sss —k k s 

which is mKdV equation. 

The biggest flaw of this model is that it is also arclength preserving which is some- 
thing that one cannot assume in the vortex-patch motion due to the filamentation 
phenomena. In fact, since the early simulations of Deem and Zabusky |DZj there is a 
strong numerical evidence (see [Arj and [Bu2]) that small perturbations of a circular 
patch can generate long and very thin fingers, and as a consequence a big increasing 
of the total length of the boundary of the patch is produced. As far as we know, 
nothing has been rigorously proved so far that can support the connection of this ge- 
ometric flow with the one of the vortex patches, (see (WD] for a discussion about this 
point and also about the connection of the travelling wave solutions of the mKdV with 
the V-states of Deem and Zabusky introduced in (DZj ). However it would be inter- 
esting to know if under some assumptions, and possibly with some extra equations 
that keep track of the evolution of the length of the curve, the mKdV can be useful 
to describe some of the phenomena that appears in the dynamics of the vortex patches. 

The rest of the paper will be devoted to the BF and its connection with NLS. In 
the next section we shall present some results on NLS equations with rough data, that 
will lead to some intuition on the existence of selfsimilar solutions for the binormal 
flow. Section §3 is devoted to the construction of some selfsimilar solutions which are 
examples of regular curves which develop a singularity in a later time. This singularity 
is either a corner or a 3d logarithmic spiral. Analogous examples but of plane curves 
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can be given for the mKdV flow, at least with some smallness condition in the cur- 
vature. We refer the reader to the paper by G. Perelman and L. Vega [PeV] for the 
corresponding results, and to the one by de la Hoz, [Ho2j for some numerical simula- 
tions concerning these examples. In the last section we shall study the stability of the 
selfsimilar solutions constructed in S3. 1-53.3. 



2 Cubic NLS with infinite energy 
2.1 Some symmetries 

Recall the cubic NLS equation 



1. l2 

iu t + u ss + -\u\ u = ^ 



y u(s,0) = u (s). 

This equation is scaling invariant: 

• If A > then u\ = Xu(Xs, X 2 t) is a solution if so is u. 
It is also galilean invariant: 

• If N € M and u is a solution then 



u N (s, t) = e- itN2+iNs u(s - 2Nt, t) (13) 



is also a solution. 
Notice that if u denotes the Fourier transform then 

un(Z) = u(Z-N). 

Hence if it belongs to the Sobolev space H k , 

Jim [\u N \ 2 (l + \H\) 2k d!i = lim [ \u(^ (1 + \^ - N\) 2k 

Jo ifk<0 
\oo if k > 0. 

Assume that we have a result that assures that if no € H k , then there exists a 
unique solution of (fT2j) up to time to = to (\\uo\\ H k) with io( r ) — > o° as r — > . Then 
if k < we will immediately obtain a global existence result from a local existence one 
just by using galilean transformations. This suggests to consider (|12p a supercritical 
equation from the point of view of galilean transformations and within the class of 
Sobolev spaces H k if k < 0. 

2.2 NLS with rough data 

In this section we want to review three results about NLS. 

First of all we want to address the issue of selfsimilar solutions with respect to 
scaling within a general framework. This was first studied by T. Cazenave and F.B. 
Weissler [CW] and then extended by F. Planchon [Pllj . |P12] to NLS equations with 
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— — — oo 

nonlinearities \u\ a u which are I? critical or supercritical, in the Besov spaces -B 2 2 • 
But recall we are interested in ID cubic NLS, which is subrcritical in I? . At this respect 
T. Cazenave, L. Vega and M.C. Vilela |CWj were able to extend F. Planchon's result 
to nonlinearities which are subcritical in L 2 although the Id cubic NLS is not included 
in their result. 

Let us consider the power like NLS 

' iu t + Au + X\u\ a u = seM. d , t6R, 

(14) 

u(s,0) = u (s), 

and focus on d = 1 and d = 2 and A = ±1. In |CVV_j the following result was proved. 
Theorem 1 Suppose 

4(d + l) 4(d + l) 

< a < 



d(d + 2) d? ' 

There exists 5 > swc/i i/tcrf i/ no 6 ^'(IR^) satisfies \\uQ\\y < 5 then there exists a 

a(d+2) , , - . 

unique solution u £ L 2 >°° (R a+i j 0/ (fllj) suc/j i/iai 

llull Q (d+2) < 25. 

L — 2 — >°°(]Rd+i) 

Here y is defined as 

y = |v? € S'OfT*) ; e itA (^ G L^ t ^ ,0Q (R <i+1 )| , 
and (|14p zs understood as the associated integral equation: 

u (t) = e itA u + i\ I e l ^ T)A \u\ a udr. 
Jo 

Notice that 

• For d = 1, — < a < 8 and we would like to obtain a = 2. 

3 
3 

• For d = 2, - < a < 3. 

' 2 

The key issue is as we see to choose a suitable space of initial data which is defined 
through the free evolution, namely 

u ey iff e itA u e L p '°°(R d+1 ), 

where L p ' 00 (R <i + 1 ) denotes the Lorentz space. These spaces are naturally defined 
by real interpolation which behaves poorly within the setting of mixed norm spaces 
L p (R, L q (R d )) and that is the reason for assuming p = q in Theorem [TJ 

The other key ingredient is to use some Strichartz estimates for the operator 



T(F) = f e^- T)A F(-,T)dT 
Jo 



(15) 



within (L P1 ,L P2 ) spaces which where unknown and involve pairs (^1,^2) which are non 
admissible in the sense of J. Ginibre and G. Velo |GVj . The possibility of the existence 
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of this type of inequalities was first observed by T. Cazenave and F.B. Weissler in 
[CW|. Based on the work of M. Keel and T. Tao [KT], D. Foschi [Fq] and M.C. Vilela 
[Vi] independently proved a general result at this respect although there are still open 
questions in this direction. 

However the inequalities for the Duhamel operator are easier to prove if one relates 
T with Bochner-Riesz operators. These operators have been widely studied in Fourier 
analysis and are given in terms of the multiplier 

™ 7 (£) = (1-I£l 2 )+, 

with (•)+ denoting the positive part. Then the Bochner-Riesz operator of order 7 is 
defined as 

BRy(fT(0 = ™ 7 (0/(£)- 
The index of interest is 7 = —1, which is closely related to the more natural 
multiplier 

#(£) = (i-|£lT\ ?eR d (16) 

understood in the sense of the principal value, and which is linked to Helmholtz equa- 
tion 

Au(x) + u(x) = f(x), x e R d . (17) 
The connection with (1151) is clear because TF is a solution of 



idtu + An = F, 

and therefore, at least formally 

~ T ~ l£l 

By scaling considerations one can always consider that supp.F(£, r) C {|£| < 1} 
and then one mimics the available proof of the corresponding inequalities for solutions 
of (fTTj) using the multiplier (fT6|) . 

The final remark we want to make about Theorem [T] is that the condition on the 
initial datum is not so easy to verify. In fact to check whether or not the condition 
is satisfied is quite a delicate question in Fourier Analysis. In the case of the free 
Schrodinger operator this question is nothing but to study LP bounds for the dual 
operator to the free propagator e %t ^UQ . In terms of this latter operator the question 
is if the inequalities 

with 

d + 2 d 
• h - = d 

p q 

• p > d, 

hold true. The answer is positive if d = 1 and it was established by C. Fefferman and 

E. Stein (see [E]). In higher dimensions d > 2 this was proved by E. Stein and P. A. 

Tomas for q < 2 (see [St]). Recently, the case d = 2 was extended by T. Tao in [Ta] 
5 

to q < — . As a consequence if [i is small enough and 

||«o|Il9,oc < " (19) 



e% n o||i>(IRd+i) < c Po|| L9 (Kd) (18) 
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for q = with — < a < 8 if d = 1 and | < a < 3 if d = 2 we can apply 

aa — 2 3 

Theorem [TJ Moreover in |CVV| it is proved that if no verifies (|19p then 

sup ||u(-,t)|| L9l00 < Cu 
t 

for a fixed constant C and with q as above. 

This was the starting point of the work by A. Griinrock [Grl| [Gr2] . In |Gr2j he 
proves the following theorem 

Theorem 2 Let 1 < r < 2 and u G L r '(R) . T/ien f/ie IKP. 

iut + n ss = A|n| 2 u n(0) = uq A = ±1, (20) 
is locally well posed and 

u(-, t) G L r '(E) (21) 

/or aZZ t. 

So in particular if 

l-o(OI < jf^p, 7 > 0, 

the conclusion holds. As it shall be explained later, the selfsimilar solutions of BF will 
lead to the case 7 = 0. 

Finally we would like to mention another work, by A. Vargas and L. Vega [VV| . 
related to (120j) and which precedes to [Gr2| . First, the following local result is proved 
in [Wj . 



Proposition 1 Assume uq G <S'(R) such that there exists T > with 

< +oo. (22) 



e itd *u Q 



L 3 ([-T,T],L 6 (R)) 

Then there exists T < T and a unique solution u of (f20|) such that 

IMIl 3 ([-T,T],L 6 (R)) < +°°- 



We say that no G 3^2 if (p2|) is satisfied uniformly in T. As a global result, the 
following theorem is proved in |VV| . 

Theorem 3 Let uo G <S'(R) that can be decomposed as uq = v o + wo with 

\\vq\\l 2 (r) - M r j Il^oll^ - f or some r <1 and M > 1. (23) 

Then there exists a unique solution u of (]20p for 

< t < T = T(M), 

that satisfies 

H u lli3([o,T],i«(M)) < c - 

Moreover 

lim T(M) = oo. 

M->oo 



10 



As a consequence if for all M > , uq can be decomposed as uq = vff + w^f 
satisfying (|23|) , then there is a unique global solution of (f20|) . 

Again the difficulty relies in giving good conditions on no such that (|22p is satisfied. 
In terms of the pointwise decay of the Fourier transform it is worse than the result 
by Griinrock, because a decay r > 1/6 is necessary. However in terms of the LP 
local integrability of uq the results are stronger. In fact it is sufficient to assume the 
following. Assume Uq 6 X p ^ which is given as the functions / such that 

i/p 

- El E 2" ro n/i r iii/ 

where T>j = {[k2 :) , (k + 1)2 J ) : k £ Z} , the family of dyadic intervals of length 2 3 

in R. Then (1221) holds p = 3, 7 = — . The above norms were firstly introduced by 

6 

J. Bourgain and naturally appear in the study of the free solutions of Schrodinger 
equation - see also |P1V| equations, (3.1) and (3.2). 

Notice that the advantage of Proposition [1] with respect to Theorem [2] is that the 
solution obtained in the former case is (R rf+1 ). This property is fundamental in 
order to integrate the Frenet equations to obtain a solution of the BF from cubic NLS. 
In any case it would be interesting to know if the methods by A. Vargas and L. Vega 
can be extended to the same range of decay as the ones obtained by A. Griinrock. 

2.3 The function as initial datum for NLS 

Assume next that uq = a5 , for a £ C. Then for all N 

uo,n(s) ■= e tNs uo = Uq. 
If uniqueness holds then we conclude from the galilean transform f|13f) that 

u(s, t) = u N (s, t) = e - itN2+iNs u(s - 2Nt, t) for all N. 
After differentiation with respect to N in the above identity we get 
= i(s- 2Nt)u(s - 2Nt,t) - 2td s u{s - 2Nt,t). 

Hence 



u(s,t) = g(t)e 
But u has to solve NLS. This gives 



2 

'77 



g(t) = — exp (±i|a| 2 hit + ia\) , a,oi £ C, (24) 
\Jt 

which has no limit as t f + . Therefore f)12[) is ill posed with no = aS (see the article 
of Kenig, Ponce and Vega |KPV2| ) . 

However notice that if we change NLS from (112ft into 



iu t + u ss + - (|n| 2 + A(t)) u = 0, 
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\a\ 2 , , s e" n 



we see that taking A(t) = , then u a (s,t) = a — -=• is a solution of 

t y/t 



1 /, |2 la' 2 



i«t + u ss + - \\u\ —j u = 0, 

u(s, 0) = a<5. 

This suggests the possibility of finding one solution of the binormal flow. The 
corresponding curve should be determined by the property that curvature and torsion 
are 

c(s,t) = ^= , T (s,t) = ^, (25) 

for some a G M + . 



3 Selfsimilar solutions for the binormal flow 

In this section we review some of the results proved in [GRV] and [GuV] concerning 
the self-similar solutions of the binormal flow. We also include in Propositions H] and 
[5] some results for curves in M 3 that as far as we know are new. These results are an 
extension to general curves of those proved in [GRV] for the particular case in which 
the curve is precisely a self-similar solution of the binormal flow. 

3.1 A V-shape as initial datum for the binormal flow 

The BF solution corresponding to (|25h should have as initial condition 




if s > 
if s < 0, 



with \A + \ = \A~\ = 1. 

Frenet equations ([8]) allow us to compute x( s ^) from (j25|) up to a rigid motion, 
obtaining a real analytic curve for t > which should have a singularity at t = . 
Notice that the binormal flow is reversible in time: if x(s, t) is a solution so is x( s > t) = 
— x(s, —t). Recall that the energy in this case is infinity, although it is locally finite as 
for the helix. 

We gather the results proven on selfsimilar solution of the binormal flow in the 
following therorem. 

Theorem 4 The self-similar solutions are parametrized by a positive number a, and 
satisfy 

\Xa(s,t) - s^4+I[ 0jOO )(s) - sA'l^^s)] < 2a\/i, 
for t > 0, where A^ are two distinct non opposite unit vectors. In particular, 

\a+s s > . , 

Xa(s,0) = \ a _ 26 
\A a s s<0. 
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Moreover, there exists a, a such that if a < a , Xa( s > 1) has no self-intersection and if 
a > a it has finitely many. Finally, the relation between a and A^ is 

9 _ na 2 

sin - = e 2 
2 

where 9 is the angle between the vectors and —A~. 

The behaviour of Xa(s, t) for t > was exhibited by the physicists M. Lakshmanan and 
M. Daniel in [LDj and a numeric study on it was done by T.F. Buttke in [Bul| , The 
rigorous description of these solutions up to t = was done by S. Guttierez, J. Rivas 
and the second author in [GRV], and will be presented in the two next sections. Let 
us mention that the link between the parameter a and the angle between the vectors 
and — A~ is an original material, and we dedicate to its proof the third subsection 
§3.3. 

Before trying to compute x is worthy to stress the point that galilean transfor- 
mations act on the set of solutions of the BF. In fact writting the flow in intrinsic 
coordinates we get -see |DaR| . |Be| 

C t = ~{ct) s - C S T 

(c ss -ct 2 \ (27) 

T ' = {—1— J, + «- 

It is straighforward to observe that for N £ R 

(cjy, t n ) = (c(s - 2Nt, t),r(s - 2Nt, t) + N) 
is a solution of (f2T|) if so is (c, r). Also 

(c n ,t n ) = (c,t) ViV iff ( CjT )=r" ±j t>0. 

3.2 The computation of the selfsimilar solutions and some 
remarks 

The simplest way to construct %(s, t) is to observe that from (125j) we can expect it to 
have a selfsimilar behaviour also with respect to scaling. Hence we settle the ansatz: 



X {s,t) = VtG{^j t>0, 



which gives from the BF equation ( xt = Xs A Xss 



-G - -G' = G' A G". (28) 
2 2 v ' 



Since G' = T , a new differentiation gives using Frenet equations 

(T A (c s n + cn s )) 
c s b + c(T A (-cT + rb)). 



S -T' = T NT" = (TA(c s n + cn s )) 



Hence (see [LD] and [Bui]), 

s 

- cn = — c s 6 + crn (29) 
2 
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and therefore 

c s = r= S -. (30) 
In particular, there exists a E M + such that c = a. Now the Frenet equations write 

V = an, 

ri = -aT +-b, ( 31) 

and we get that G' = T is a real analytic function. Moreover 

T(M) = r(A). (32) 

Recall that BF is invariant under rotations so we can assume without loss of gen- 
erality that 

t\ It 

n ](0,0)=I 3 A3 and from (1321) ( n |(0,t)=I 3 A3- 

We have to determine x(0> t) . But from 

Xt = cb 

we conclude 

Xt(0,t) = -j= (0,0,1). 

As a conclusion x(0, t) = 2ay/t (0, 0, 1) . We could have got it also from ([28 
Therefore 



and G is determined up to a rotation by the solution of (|3ip with 

G(0) = 2a (0,0,1) 
T(0) = (1,0,0) 
n(0) = (0,1,0). 

It is easy to do the numerical analysis of G . We start with (|28p : 

-G - -G' = G'A G" 
2 2 

Recall the identity: 

F A (G A H) = (F ■ H) ■ G - (F ■ G) ■ H. (33) 

Then 



G' A [G' A G") = (C ■ G") G' - (G' ■ G') G" 
But on the other hand from 

~2 
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S -T' = T A T" 



we get 

and \T\ 2 = 1. 

Hence from ([28]), (G2 



-f T' T = = --G' G" 
2 2 



f -GAG' - G" (34) 
1 G(0) = (1,0,0) G(0)-G(0) = 0, 

which can be easily computed numerically. We refer the reader to the web page 
http: / / ehu.es /luisvega/ vortex/files / esquina.mov to see the corresponding flow of curves. 
In these pictures is clearly seen that the profile G tends to two straight lines at infinity, 
and as a consequence that 

, . f A + s s > 
40= . 35 
\As s < 0. 

In order to prove it we will have to work a little bit. We follow the arguments of S. 
Gutierrez, J. Rivas and L. Vega [GRV]. Recall that from ([28]): 

~G — - G' = T AT' = T A an = ab. 
2 2 



Hence, for s > 



It follows that 



G(s)\'_sG'-G_ 2a 



s 2 s 2 



G(s) 



f' s 2n 

/ -^b(s')ds', (36) 
S Ji s u 

and therefore, since b is unitary, we conclude that there exists 

l im ^£l = A+. (37) 

s— >oo s 

G(s) 

Analagously we obtain that there exists A~ such that lim = A~ 

Moreover, 



and 











s— >— oo s 


A+- 


G(s) _ G(s) 
s s 


oo 
s 


/•oo 

-2a / 

J s 


b{8,) ds> 


G(s) 


POO 

= sA + + 2as J 

J s 


s' 2 


s > 


G(s) 


= sA~ - 2as [ 

J — oo 


s' 2 


s < 0. 



(38) 
(39) 



We have the following proposition. 
Proposition 2 For all s 

\ X (s,t) - sA + I [0iOo) (s) - sA-I { ^ Q] (s) \ <2aVi t > 0, 
where I denotes the characteristic function. 
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Proof If s = it is trivial because %(0, t) = 2ayi. We will assume that s > 0, the 
other case been proved in a similar way. Then 

(s, t) = VtG ( 4=1 = sA+ + 2as /°° ds', 



M s 
6(s' 



M *' 2 



■ (is' 



X 

and then 

\x(s,t) - sA + \ < 2as 

The proof follows using that b is unitary. 
Remark 1 It can be also proved that (see IGRV)j ) 

, ™, \ M + s > 0, 
UmT«, = 

t-+o v 7 s < 0. 

Our next purpose is to prove that A + ^ A~ . Recall that we are dealing with a 
one parameter family of solutions with a as the free parameter. We will make explicit 
the dependence on a in the following proposition in order to clarify the exposition. 

Proposition 3 (i) A+ j= -A~ Va. 

(ii) A+ = A~ <=► a = 0. 

(Hi) 3 a, a such that if a < a , G a (s) has no selfinter section and if a > a it has 
finitely many. 

Proof Recall that G a - s G' a = 2a b. Also G" = T = an. Hence 

G^ = an' = a^-aT + S -b} 

= -a^T+ S -(G a -sG' a ). 



Then 



s 



G:+^+--^G' a - S -G a = (40) 
G o (0) = 2a (0,0,1) 

G' a (0) = (1,0,0) (41) 
G'M = a (0,1,0) 



Also, 



x a (0) \ / 

G a {s) = (x a {s),y a {s),z a (s)) with | x' a (0) = 1 
• In ([4*0]) if G a (s) is a solution then G a (— s) is a solution. 
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We conclude that if a ^ then 

Xa(s) = -X a (-S), 

Va(s) = y a {-s) : (42) 

Za{s) = Z a (s), 

and therefore x a , y a , z a are three linearly independent solutions of the scalar ODE 



w" + [a 1 + — j w - - w = (43) 

which are tempered distributions (G a ,G' a ,G'^,G'l' grow polynomialy). Applying the 
Fourier transform we get: 

£Q" + 3Q' + 4£ 2 w -4£a 2 £ = 0. (44) 

This is a singular regular O.D.E. with r = and r = — 2 the solutions of the indicial 
equation 

r(r — 1) + 3r = 0. 

Therefore the regular solutions have to be even, and the odd solutions have to 
behave as uiq 5' at the origin. But we know that x a is odd and non trivial. Therefore 
w a which amounts to say that A^i ^ 0, where we use the vectorial notation 
At = (A± v Aj 2 ,A± 3 ). From g2D we get 



a,l' 



4J"i = A 

■A£,3 = ~A a ,3- 

Hence A+ + A~ ^ 0. 

Assume now that A+ = A~. Then A+ 2 = A~ 2 = A~ 3 = A+ 3 = 0. But then 
the even functions y a (£)> ^a(C) should be regular at the origin and linearly dependent. 
Hence a = and = (1, 0, 0) =A~. 

Only (iii) remains to be proved. Recall that G a = sT + 2ab. Hence |G a | 2 = s 2 + a 2 . 
Therefore if G a {s\) = G a (s2) for different s\ and S2, then s\ = s\ and si = — S?,- But 
from (|42p this implies that a;(si) = = x(s2). Recall that 

x "' + {f + t) x ' a ~\ Xa = °' 

Xa(0) = 0, 

<(o) = i, 
<(0) = o. 

X { S } 

Hence x a (s) = s + O (a 2 s 3 ) and on the other hand the map a i — )• 

lim = ^+ 

s— >oo s 1 

is continuous from the proof of Proposition [2J Hence if a is small enough there is no 

/s 3 

selfintersection. On the other hand it is easy to see that x' a (s) = cos(as) + O I — 



and therefore 

^o(s) = - sin(as) + O ( — ) < s < 1. 
a V " 



a 
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Hence for a large there are selfintersections. 

Notice that if x( s ^) is a solution — x(s, —t) is also a solution. Hence 

Xt = Xs A Xss 
X (s,-1) = -G a (s). 

is a solution that develops a singularity at time t = although G a is real analytic. 
A few remarks are in order. 

Remark 2 The kinetic energy (J c 2 (s, 1) ds) is infinity but locally finite as in the case 
of the helix. 

Remark 3 The solution Xa (s,t) = \ftG a (^y^j ^ as as t° r sion — which is selfsimilar 
of the type 

T X (s,t) = \ T(\s,\ 2 t). 

This scaling leaves invariant the following conserved quantity 

T(s ,0M» = /r(MM, 
Notice that in our situation the above invariant has to be understood in the sense 



f 

lim / t(s, t) ds = 0. 
R-+00 J_ R 



Remark 4 It is interesting to analyze the type of singularity we get from the point of 
view of Euler equations and vortex filaments. In particular let us recall the condition 
obtained by P. Constantin, C. Fefferman and A. Majda \CFM$ which is measured in 
terms of \\^ x £,\\l°° , where £ stands for the unit direction of the vorticity: 

Their condition is that 

I HVa.fHi.dt = +oo (46) 
Jo 

if there is a blow up at time T. In our setting f plays the role ofT and therefore \V X £\ 
is the curvature. Hence ([46 p becomes 

a 2 

— dt = oo. (47) 

Remark 5 Recall that in the case of the free Schrodinger equation we have a similar 
situation. Taking 

e 4* 

u(s,t) = —j= - 1< t < 0, 



we have a solution of 

( iu t + u ss = 0, 



.2 



u(s,-l) = e J 4, 



such that u(s,0) = 2y / 7re*4 S. 
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Remark 6 Equation (jlOp of the tangent vector suggests many possible generalizations 
by considering other targets besides the sphere. Let us then introduce the notation 



10 
1 
0-1 



u A v. 



Therefore instead of (|10|) we consider 



T t = T A- T, 



SS 1 



(48) 



with T a map from M? onto the hyperbolic plane M 2 . Analogously, since T = \s, we 
can obtain the equation 



Similar calculations and results as in JGRVf were done by de la Hoz JHol\/ for self similar 
solutions of (|49p in the hyperbolic setting. The extra- difficulty is that there is a-priori 
no control on the size of the euclidean length of the generalized binormal vectors. 

Remark 7 Let us finally remark that numerical computations on the selfsimilar curves 
of BF were done by F. de la Hoz, C. Garcia- Cervera and L. Vega in \HG-CY$ . In that 
paper the authors point out the similarity of these curves with those that appear in the 
experiments of the flow around a delta wing. 

3.3 Proof of the identity A a \ = 

Let us try to answer a more delicate question. We have been able to construct non 
trivial solutions of BF such that the initial condition has a V-shape determined by 
two unit vectors and A~ . The question we propose is if given any pair of vectors 
(A + ,A~) there is a such that (A + ,A~) = (A£,A~). We have seen in Proposition [3] 
(ii), that the particular case A + = —A~ cannot be obtained using the family of curves 
{G a }, a£l. In this section we shall prove that except for this case the answer to our 
question is positive. 

The proof requires to integrate the Frenet system of equations in a different way. 
We start defining the complex vector 



Xt = Xs A_ x. 



(49) 



N = n + ibe i ^< s ') ds ' 



and 



u = ce J o 



.* Jo ^O')^' 



Then we rewrite the Frenet system as 



N 3 = 




T s = 



We denote (vi,V2,Vs) the components of a M 3 vector v. Then 



ni + ibA 2 + T? = n 2 + b 2 




= 1. 
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Next we introduce 



so that (fj satisfies the Riccati equation 

/ 1_ 2 1 

Now we consider r\j such that 

. U . . U lol 1 , 2 1N 

= c^i + lT ~^ = 2 U7] J + 2 U= 2 U W + ' 

Hence ^ 

Vj + irr/j = -c (l + Vj) ■ 

The next step is to define Oj as Oj = exp[— c/2 rjj], i.e.: 

2 «' 

Then 

c' # 2 0" 2 0' 2 

n'. = 9 i- -I J — 

Vj c 2 0j cOj c 0) 

Hence 

2 c 2 % c%~ iT c9 j + 2 C ' 
We have arrived at the remarkable fact that Oj solves 

( -^ + iT )o' + ^0 i = o. 



5 

e 2 



Multiplying by we get that 



(Ql \ z Q2 
-j-j constant = E(0). 

Also notice that 

Wj\ 2 + \Tj\ 2 = l 
is equivalent to (rij — ibj)(rij + ibj) = (l — Tj){l + Tj). Hence 

l9 rij+ibj rij—ibn l + 

and 

_ |^| 2 -l 
20 



. p 2(9' 
But (pj = e Jo T r]j and rjj = — -p Then 



' " | % | 2 + 1 fl^ + l^-l 2 ^ 2£(0) 



Therefore we have proved the following proposition valid for any curve. 
Proposition 4 Any of the components Tj of the tangent vector T can be written as 



\0i\ 2 
Tj = 1 ! 1 



3 * 2£(0) ' 
with 8j solution of [50j) and E(0) as in [51\) . 

Assume for simplicity 

yn J (0) = w 

Then 

ri(o) = i =>- 0i(o) = o. 

c(0) 



Taking 6>i(0) = ^ we get 



T\ = 1 — |^!| 2 ; £(0) = ~. (52) 



Similarly 

T 2 (0) = =>- %(0) = 1 , 2 (O) = 1 , # 2 (0) 



c(0) 



Then 



Finally 



gives 



T 2 = 1-|0 2 | 2 ; £?(0) = i. (53) 

T 3 (0) = , 773(0) = i , 9 3 (0) = z , 0> 2 (0) = ^ 

r 3 = i-|e 3 | 2 ; E(0) = ±. (54) 
As for n and b we get, dropping the j -index, 

• cn = T s = (l - |fl| 2 ) s = -2Re6>#' 

2 / 

• T ss = -2\9'\ 2 - 2Re69" = -2\e'\ 2 + — \6\ 2 - 2-Re9e' + 2rlm00 / 

2 c 

From T ss = (cn) s we get 

-b = lm9e'. 
2 

Hence we have proved the following proposition valid for any curve. 
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Proposition 5 Any of the components nj and bj of the normal and binormal vectors 
n and b can be written as 



c(nj — ibj) = —28j0j , j = 1,2,3, 
with 9j solution of l50j) verifying {oT^) -(54\). 

In our particular case equation ([50]) reduces to 



(55) 



2 4 



(56) 



Our final purpose is to find a relation between 7, the angle determined by and 
—Aq, and the parameter a. 

It is easy to find two linearly independent solutions by just computing the Fourier 
transform in 



hi?) 



ft 00 



-^ 1Ogl V + oc)(0 



e i(sz+e)A. 



Recall that A+ = (A aj i, A a>2 , A a>3 ) , A a = (Ai,i> -A a>2 , -A a)3 ) so that 

2 1 
cos 7 = 1 - 2A al ; sin - = A a>1 . 

Then it suffices to compute A a 1 = lim Ti(s). 
But 



r x (s) = 1 - |^i| 2 = 1 - louA + 021ft 



with 



Hence 



= «2,1 



2^2ft(0)' 



1-Tx 



IW IAW+AWI ' 



But notice that 



and 



lim (ft+ft)( s ) 



d£ 



e 2 



i%- log |£| 



d£ 



ft(0) 

After some calculations we get 



2 fOO 



A a ,i = e 2 . 

Finally notice that the dispersive relation of ft and ft is 

<K£)=£ 2 -ylog|£l, 



(57) 
(58) 



(59) 



(60) 



which indicates that the non-linearity only makes a logarithmic correction to the free 
evolution. 
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3.4 Spirals 



In this section we will review some of the results obtained by S. Gutierrez and L. Vega 
in |GuVj . So far we have dealt with the I.V.P. 

iu t + u ss + ^{\u\ 2 + A{t)) = 

(61) 

u(s,0) = a 5. 

In this section we will be interested in 

u(s, 0) = ai p.v. -. (62) 
s 

On the other hand concerning the BF we have characterized all the possible solu- 
tions with respect to the scaling that preserve arc length. They were determined up to 
rotations by a free parameter, namely the curvature must be —p and the torsion — . 

In order to capture the solutions of the BF which come from the initial condition 



(|62h we have to change the ansatz. Assume A is a real antisymmetric matrix. Then 
consider solutions of the BF 

Xt = Xs A Xss (63) 

of the type 

x ( 8 ,t) =e A / 21ost ViG(-^\ t>0. (64) 

Hence G has to be a solution of (e-^/ 21og * is a rotation) 

(I + A)G — sG' = 2G' A G" , \G'\ 2 = 1 . (C ■ G" = 0) . (65) 
Multiplying on the left by G' we obtain 

G" = ^{I + A)G A G', (66) 
which together with the conditions 

(I + A)G(0) -G'(0) = 0, (67) 

|G'(0)| 2 = 1, (68) 



is equivalent to (|65p . It is quite striaghforward to prove that (|66p has global existence. 
Also it can be easily computed numerically. Again in this case the relevant feature is 
the behaviour of G at ± infinity. This follows from the identity 

e -Aio gs G(sy\ _ {-A + I)G-G' _ 2 cb 
s J s 2 s 2 

Hence if we prove that the curvature is bounded we can proceed as in the case of 
the kink solutions. This is a consequence of the following elemental lemma. Consider 
without loss of generality 











/' 


















A= (i fi G 
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Lemma 1 If G solves (|65p then 

\T'\ 2 (s) = -fiT 3 (s) - v with (69) 

v = -^T 3 (0)-±\(I + A)G(0)\ 2 . (70) 

Proof Differentiating in (|65p we get 

AT = sT' + 2T A T'. 

Then 

AT A T = -s T A T' - 2 (T • T') T + 2 T" 

(iT A T) • T' = 2T" • T = 4- \ T '\ 2 - 

as 

But (AT A T) ■ T' = —fiT^ which concludes the proof. 

/ AT AT A T \ 
It turns out that if a ^ the frame I T, , — — , — - — — ) is the right one to study 

r V A r|y 

this family of solutions in the spirit of what we have seen for the kink solutions. This 
is a consequence of the identities: 



ds \ tJ J \\AT A T| 2 |AT| 

Let us consider (y, /i) given by 



cfc 2 dlT'l 2 



,,.~-S(r-'-y (73) 



The key point then is that (y, h) solves the system 
' y' = sh + g(\T'\ 2 ), g(\Tf) = 2E(0) 

h' = - S -rV, 



(3\T f \ 2 + v)(\T'\ 2 + v) 
-yu j i j> <ai j i i = m\y> 

s 



(74) 



with ^ given in (|70l) . 

The relation with NLS of the above comments is also clear. The ansatz ()64j) implies 
that c and r have to be selfsimilar i.e. 

Therefore we have to look for solutions of © with A(t) = — of the type —= u ( — 

t y/t \ v t 

This gives us ODE 

u" - l -(u + su') + ^(\u\ 2 + v) =0. 
24 



Then we consider u(s) = f(s)e ls2 ^ 4: with 



r+^r+^(i/i 2 +^) = o. (75) 

The fundamental property of (|75p is that it has a natural energy 

E(s) = \f\ 2 + 1 -{\f\ 2 + u) 2 = E(0). (76) 

Finally the way to analyze the solutions of (f75|) is by defining y and h as 

l+ih = Jf. (77) 

It turns out that (y, h) solves the system (J7H). The results can be found in [GuV| . 
With respect to the I.V.P. with data (f62|) we have the following result. 

Proposition 6 Let v > and consider the I. V.P. 



U /, , 2 V 



iu t +u ss + - (Jtip + ~) =0 t>0 sel, 

u(x,0) = a\p.v.- a\ G C\{0}. 



(78) 



T/ien either there is no weak solution u for the I.V.P. (|78p m £/ie c/ass 

i^V t 3 / 2 |n| 2 uGL°°([0,+oo),5( 



wzi/i lim u(s,t) = a\p.v. ( — ) in S'(M.) or there is more than one. 
40 V s 1 



4 Stability 

In this section we shall present the known results concerning the stability of the BF 
selfsimilar solutions Xa, which show that under some appropriate small regular per- 
turbations the corner still remains at time zero. To this purpose we first study the 

perturbations of the filament function u a (s,t) = a^E-, solutions of 

1 (^ ,2 V?\ 

iu t + u ss + - \ \u\ — I u 

u(s,0) 

The formation of a corner for the flow of curves Xa corresponds to the limit a 5 oiu a as 
time goes to zero. The informations about the perturbations of u a will be gathered in 
the first subsection. Then in the last subsection we shall deal with the construction of 
new BF solutions. We shall define from a perturbation of u a a curvature and a torsion 
function. Then we shall construct a flow of curves, solution of the BF with these 
curvature and torsion function. Finally, we shall prove that the new flow of curves is 
close to the selfsimilar one, and in particular a singularity is still generated at time 
zero. 



(79) 

= a 5. 
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4.1 Perturbations of NLS solutions with ^-function as ini- 
tial datum 

In an attempt to understand the above equation (fT9|) . we considered in |BV1] the 
equation 

iu t + d ss u±\u\ a u = 0, (80) 
with an initial data of Dirac type, and its solution 

u(s,t) = u a)±a (s,t) = u a (s,t) e ±iAa - a{t \ 

with 



\a\ a • / r, 



lap log t si a 



For a > 2, u a ^± a does not have a limit as t goes to zero, and this gives an instability 
result, as explained in the previous subsection 2.3. 

A first natural question is to find out if u a ^± a is a stable solution in the case a < 2. 
Denoting 

r l (s,t) = e^ A ^ t \u-u ai±a )(s,t), 
and using the fact that u a is a solution of the linear equation, we have 

f ir\ t + d ss r\ ± (\r\ + u a \ a - \u a \ a ) (r? + u a ) = 0, 
\ r/(s,0) = u (s). 

If a < 2 then the term \u a \ a is integrable at i = 0, and we can perform a fixed 
point argument. We have obtained the following positive answer, valid also in higher 
dimensions, for a < § (a similar result is presented in |BV1| for Sobolev spaces, with 
a more technical proof). 

Theorem 5 Let a < 2 and uq G L 2 (M). There exists a time T = T(a, \\uo\\l 2 ) an d 
a unique solution u of equation (180p with initial data u$ + a 5, such that 



u - u a , ±Q e L 4 ([0, T), L°°(R)) n C ([0, T), Z/ 

A case of rough perturbation of o5q has been studied by Kita. In [Ki] he has 
described the solution generated by the sum of two or three Dirac masses. 
Finally, we shall do a parallel with the corresponding parabolic problem, 

ut — Au ± \u\ a u = 0, 
u(s, 0) = a 5, 

which has been intensively studied. F.B. Weissler has shown in |We| that in the fo- 
cusing case, for a < 4, there is no unicity of the solution. In the defocusing case, if 
only positive solutions are considered, H. Brezis and A. Friedman have proven in [BFJ 



the existence of a unique solution for a < |, and non-existence for a > 4. Therefore 
for the heat equation there is an important difference between the focusing and the 
defocusing case, even for a < |. Theorem [5] shows that this is not the case for the 
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Schrodinger equation. 



The second natural question is to find out if u a is a stable solution of the cubic 1-D 
Schrodinger equation modified by the complex phase e ±m log *. 

Having as a purpose the study of perturbations of the selfsimilar flow Xa, we shall 
first consider perturbations at positive times of the filament function 

u a (s,t) 



>/i' 

and we shall treat more generally the equation 

iu t + u ss ± \\u\ — — I u = 

around its particular solution u a . By defining v the pseudo-conformal transformation 
of u, 

e iS & ( s 1\ 

" (M) =vT <81) 

this is equivalent to study the large time behavior of perturbations of the constant 
solution a of 

iv t + v ss ±-(\v\ 2 -a 2 )v = 0. (82) 
Let us notice that this equation has a natural energy 

E{t) = \ [ M*)| 2 ^^ f(\v(t)\ 2 -a 2 ) 2 ds. 



with the law 



d t E{t)T^ J(\v\ 2 -a 2 ) 2 ds = 0. 



In the defocusing case, this has allowed us to obtain in |BV1| a good control in time 
of the norm \\v(t) — a||^2 and to conclude global existence. 

As the coefficients of the linear terms are of the type j, and since the perturbation 
is done around a constant function, there will be connections with the long range 
effects of 1-D cubic Schrodinger equation (T. Ozawa [Ozj . R. Carles |Caj . N. Hayashi 
and P. Naumkin |HNJ) and with the scattering around the constant solution 1 of 2-D 
Gross-Pitaevskii equation (S. Gustafson, K. Nakanishi and T.-P. Tsai [GNTJ). 

We shall search for wave operators for the perturbation w(t) = v(t) — a, meaning 
that we shall search for solutions of 

iwt + w ss db - (jw + a\ 2 — a 2 ) (w + a) = 0, 

behaving at large times as a given profile w\(s,t). This can be done by a fixed point 
argument, in a suitable space defined around w\, for the application 

Aw = WiTi f°° e i { t- T )dl f (\w + a\ 2 -a 2 )(w + a) _ (K + a\ 2 - a 2 ){ Wl + a) \ ^ 



t* / e *->« ( (K + "l -a )(«>! + a) _ {ldr + dl)wi , dT 



2\ 
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For showing that w behaves like a free evolution e ttds u + , one should control the 

2 

Duhamel terms corresponding to the worse source terms, the linear ones ^j-w\ and 

2 

Qj-wx- These Duhamel terms are 



e i{t-T)8* jrdl u _ 



dr 

T 



/oo 



dr 



We notice that the first integral in not controllable. To avoid this term, we shall look 
for w close to the long range ansatz 

Wl (s,t) = e ±ia2lost e iW ° u + (s). 

By using inhomogeneus Strichartz inequalities we manage to control the remaining 
linear term 



l±ia 2 



dr 



Li 



-l±ia 2 



<Ca 2 



Co 2 \\u + 1 



Lp' ((t,oo),Li') 



Li' 



dr 



1 



Li 



Lp' (t,oo) 



Ca \\U-\- H^q' — . 

tp 



The 1-admissible couples (p, q) are between (oo, 2) and (4, oo), therefore if u + G L 1 nL 2 , 
the best decay we can obtain this way is t~ 4. We denote 

Vl (s,t) = a + wi(s,t) = a + e ±ia2logt e it9 ° u + (s). 

We get our first result of [BV2j . in mixed norm spaces L p L q . 

Theorem 6 Let to > 0. There exists a constant ao > such that for all a < ao, 
and for all u+ small in L 1 n L 2 with respect to ao and to, equation (\82h has a unique 
solution v with 

v - Vl 6 C ( [t , oo) , L 2 (M) ) n L 4 ( [t , oo) , L°° (R) ) , 
satisfying, as t goes to infinity, 

\\v{t) - Vi(t)\\ L 2 + ||v - Ul||L4(( t ,oo),£'») = 0(t~%). 

However, the t~i decay, optimal for the methods used in mixed norm spaces, is not 
enough to construct BF solutions. We shall then exploit more the oscillations of the 
linear remaining term, by using the Fourier transform and strengthening the conditions 
on u+. Writting 



i(t-2r)d s 2 I r -iTdj "+ 
T-l±ia 2 



dr 



and doing an integration by parts in r we get 



Using Plancherel formula we obtain the control 

-ird 2 



~-l±ia 2 

e -i{t-2r)e 
~2±ia 2 



e ix Z t£(-f)d£dr, 



Jst, 



2ie 



d^dr. 



J(t-T)d 2 



-l±ia 2 



dr 



<"1 



L 2 



£P(-0 



L 2 



The derivatives V fc of this term have a 1/i decay, with constants depending on the 
norm Q f obtain in [BV2| a new result, with a better decay. 
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Theorem 7 Let to > and m € N* . There exists a constant uq > suc/i i/iai for all 
a < oq and for all u+ small in H~ 2 n H m n W™ 1 ' 1 with respect to ao and io> equation 
has a unique solution 



v- Vl €C([t ,oo),H m (R)), 
satisfying, as t goes to infinity, for all entire number < k < m, 

\\{v-v 1 ){t)\\ L ,=0{r^) , \\& a (v-v 1 ){t)\\»=0(t- x ) 



(83) 



The presence of two different rates of decay is now due to the quadratic terms 
^|f«i| 2 and jw\. The less oscillating is the first one, which gives the rate of decay t~2 . 
However, its derivatives can be better upper-bounded. We shall sketch in the following 

I |2 

the way we treat the Duhamel term coming from the worse quadratic term , 1 V . We 
shall use again the Fourier oscillations, 

|2 X 



d k I e i(t-r)d] 




lei' 



<t-r)e 



\t\ k e is ^{r ] ) : u+{r ] -Qdr ] d£,dT, 



and an integration by parts in r 

-it{e+2£,{r,-H)) 



t 



\ee 



k is£ 



2i£(v - 



dr\ d£ 



+ 



K|fc c M u+(r])u + (r] - Q 



drj d£ dr. 



Hence we have obtained that this quadratic Duhamel term can be controlled in H l by 



Ca 



U+ -k 



L- 



Ca 
T 



L- 



Ca u , . , 
< —\\ u +\\hA\ u +\\h-^- 



We notice that this argument can be repeated for H k , with k > 1, but not for k = 0. 

At the level of the initial equation before the pseudo-conformal transformation, we 
obtain that u(t) is close to 



ui(s,t) = u a + e 



The condition u + € H~ 2 is translated on u + by a canceling condition at the origin, 
which could be related to the blow-up result of J. Bourgain and W. Wang [BWj. We 
notice the absence of a limit in I? for u(t) — u a at time zero. However this does not 
exclude the possibility to construct from u a binormal flow with a trace at zero. This 
is because the curvature and the torsion defined from u corresponds to the second 
derivatives of the curves. 
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4.2 On the stability of selfsimilar solutions of the binor- 
mal flow 

In this last subsection we sketch how from a solution of Theorem [7] we construct a flow 
of curves Xi solution of the BF, and get information about its limit as t goes to zero. 
We prove that the new curve \ is close to the selfsimilar one Xa, and then obtain the 
persistence of the formation of a singularity at time zero. 

Let v be a solution of Theorem [7] with m = 3, and u its pseudo-conformal transform 
(|8ip . Since v is regular enough and does not vanish, we can define two real functions 
c and (j) such that 

u(s,t) = c(s,t)e i4 '^ t \ (84) 

We define r(s,t) := (p s (s,t), so 



u( S ,t) = c(s,t)e t fo T{s ' t)ds+ ^ (} ' t) . 

c(s,t)e iJ ° T(s,t)ds 



Then, the function 

u(s, t) 

is a filament function, solution of ([5]) with A(t) replaced by + <pt(0,t). In partic- 
ular, (c, r) satisfy the intrinsic equations (|27[) . and a computation yields the extra- 
information at s = 



T +^(0,t) 



CT 



(0,i)+c 2 (0,t). 



(85) 



In view of the pseudo-conformal transform, the curvature is 



c(s,t) 



, \ i 1 
= 7^ 









i) 







and the torsion is 



r(s,t) 



u s{s ,t) _^v(l, 1 J )+d s v(l,\) 



As v £ H 3 (R) C 
By definition, 



, the curvature and torsion function will be regular enough. 



v(s, t) = a + e 



±ia 2 log tjtdi 



u + (x) + (v - vi)(s,t). 



The decay of Theorem [7] combined with the explicit formula of the free Schrodinger 
evolution allows us to obtain that (c, r) is close to (c a ,r a ), with strong rate of decay 
as t goes to zero, and constants depending on the norms of u + . 

Starting from (c, r) we construct the binormal flow as follows. We fix an initial 
data (T, n, 6)(0, to) at time to = 1/to, to be specified later. For to > t > we define 
(T,n,b)(0,t) by imposing that the derivatives in time satisfy the same system as the 
tangent, normal and binormal vectors of a binormal flow of curves satisfy in general. 
Then, we define (T,n,b)(s,t) from (T,n,b)(0,t) by integrating the Frenet system for 
fixed time t. We therefore impose 



T 

n 
b 



(s,t) 




(0,*o)- 



CT 



-CT 





\ 



n 
b 



{0,t')dt' 
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— r 


o y 







+ 

Proceeding this way, T shall satisfy the equation 

T t = TAT ss , 

of the tangent vector of a general binormal flow. 

Once T is constructed, for a given x(0, £o) we define 

X (s,t) = x(0,to)- [\cb)(0,t')dt' + I T(s,t)ds. 
Jt Jo 

By using Frenet system, 

T t = T A T ss = T A (c n) s = T A (c s ra + c r 6) = -crn + c s 6, 

and it follows that x satisfy the binormal flow. 

Finally, t) is obtained this way for all times when (c, r) are regular, that is 
to > t > 0. Knowing that c is close to c a = we get 



Ix(Mi) - X(s,i2)| 



*2 



c(s, t) b(s, t) ds 



< 



* 2 Ca , 
— =d£ — >■ 0. 

1 tl,*2^0 



This implies the existence of a limit Xo( s ) at t = 0, and we obtain similarly that for 
all s £ (— oo, oo), 

|xM)-Xo(s)| < CaVi. 

Next, for proving that x( s )0) is close to Xa(s,0), we first show that the tangent 
T(s,t) is close to T a (s,t). As a first step, we prove that T(0, t) is close to T a (0,t), and 
then we prove that T(s, t) is close to T a (s, t). For the second step we need to know that 
(T, n, b)(0, t) is close to (T a ,n a , b a )(0, t). In the following we shall sketch the estimates 
at (0,t), which is the less technical step, but nevertheless a key point of the proof. 

From the above definition of (T, n, 6), 



n 
b 



(0,t) 



/ o 

CT 



-cr 




(^) 



c 






Now we use the extra-information at x = (|85p . 

.2 



cr 



Therefore we obtain 



/ T N 
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-cr c s 
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cr 
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In order to get rid of the term (j) t (t,0), we perform the change of functions h + ib = 

e l 2(n + ib), and we impose as an initial condition (T, n, b)(0, to) = (T a ,n a ,b a )(0,to). 
For having (T, h, 6)(0, t) close to (T a ,n a ,b a )(0,t), it is then enough to show that the 
matrix 



/ \ CT \ \°s\ \ 
\ct\ 



V 



2 





(t,0) 



J 



2 

is integrable in time near t = 0. This is true, in view of the strong rate of decay 
obtained for (c, r) — (c a ,T a ). Also, we get a good decay for |e*2 — l| 5 so (T, h,b)(0,t) is 
close to (T, n, 6)(0, i). In conclusion, we have (T,n,b)(0,t) close to (T a , n a , b a )(0, t). 

We obtain finally that for e > 0, and u + small enough, 

|(r-r a )(M)|<€. 

In particular, for s > 0, 

ps ps ps 

X(s,t)- X (0,t)= T(s,t)ds = A+s+ (T-T a )(s,t)ds + T a {s,t) - A+ ds, 
Jo Jo Jo 

The third term tends to zero as t goes to zero, uniformly in s |GRVj . By letting t go 

to zero, we obtain the desired information 

|Xo(s)-Xo(0)- a/4+1 <es. 

More precisely, we have obtained in [BV2] the following result. 

Theorem 8 Let e > 0, to > 0, < a < do, where ciq is the constant in Theorem^ 
Let u+ small in H~ 2 n H 3 D W 3,1 with s 2 Uj r small in H 1 with respect to e, a and to, 
and let v be the corresponding solution obtained in Theorem^ By using the Hasimoto 
transform, we construct from v a family of curves x( s ^) which solves the binormal 
flow for to > t > 0, and such that there exists a unique xo satisfying 

|x(*.*)-Xo(«)| < CaVi 

uniformly in s E (—00,00). 
Moreover, 

|XoO) - Xo(0) - sA+I [0jOo) (s) - sA'I^^^is^ < e|s|, 
where A^ are two distinct non opposite unit vectors, such that the angle 7 between A+ 
and —A~ is determined by the relation 

7 7ra 2 

sin — = e 2 . 
2 

Since 

Xa(s,0) = s J 4+I [0jOO )(s) + s^~I ( _ OO)0 ](s), 

we have that xo( s ) nes hi the e-cone around x o (s,0). Therefore we have obtained a 
large family of curves, solutions of the binormal flow, which are close to the selfsimilar 
solution Xa- In particular a corner is still formed for x( s >t) at t = and at s = 0. 
The angle of this corner can be made as close as desired to the one in between A^ and 
—A~ by taking e small enough. This shows that the formation of a singularity in finite 
time of Xa is not an isolated phenomenon. In a work in progress we study, instead of 
the existence of wave operators in Theorem [71 the asymptotic completeness question, 
which could lead to a better description of the allowed perturbations. 



32 



Alekseenko, S.V., Kuibin, P.A., Okulov, V.L. Theory of concentrated vortices: 
An Introduction, Springer Verlag, 2007, ISBN: 978-3-540-73375-1, 494p. 

S. V. Alekseenko, P.A. Kuibin, V.L. Okulov, S.I. Shtork, Helical vortices in 
swirl flow, J. Fluid Mech. 382, (1999) 195-243. 

H. Aref, Integrable, chaotic, and turbulent vortex motion in two dimensional 
flows, Ann. Rev. Fluid Mech. 15 (1983) 345-389. 

V. Banica, L. Vega, On the Dirac delta as initial condition for nonlinear 
Schrddinger equations, Ann. I. H. Poincare, An. Non Lin. 25 (2008), 697- 
711. 

V. Banica, L. Vega, On the stability of a singular vortex dynamics, Comm. 
Math. Phys., to appear. 

R. Betchov, On the curvature and torsion of an isolated vortex filament, J. 
Fluid Mech. 22 (1965) 471. 

J. Bourgain, On the restriction and multiplier problems in R 3 . Geometric 
aspects of functional analysis (1989-90) , 179-191, Lecture Notes in Math., 
1469, Springer, Berlin, 1991. 

J. Bourgain, W. Wang, Construction of blowup solutions for the nonlin- 
ear Schrddinger equation with critical nonlinearity. Dedicated to Ennio De 
Giorgi. Ann. Scuola Norm. Sup. Pisa CI. Sci. (4) 25 (1997), 197-215. 

H. Brezis, A. Friedman, Nonlinear parabolic equations involving measures as 
initial conditions, J. Math. Pures Appl. IX, Ser 62 (1983), 73-97. 

T.F. Buttke, A numerical study of superfluid turbulence in the Self-induction 
Approximation, J. of Comp. Physics 76 (1988) 301-326. 

T.F. Buttke, A fast adaptive vortex method for patches of constant vorticity 
in two dimensions, J. Comput. Phys. 89 (1990), 161-186. 

R. Carles, Geometric Optics and Long Range Scattering for One-Dimensional 
Nonlinear Schrddinger Equations, Comm. Math. Phys. 220 (2001), 41-67. 

T. Cazenave, L. Vega, M. C. Vilela, A note on the nonlinear Schrddinger 
equation in weak L p spaces, Comm. Contemp. Math. 31 (2001) 153-162. 

T. Cazenave, F.B. Weissler, Asymptotically self-similar global solutions of the 
nonlinear Schrddinger and heat equations, Math. Z. 228 (1998) 83-120. 

P. Constantin, C. Fefferman, A. Majda, Geometric constrains on potentially 
singular solutions for 3-D Euler equations, Comm. P.D.E. 213&4 (1996) 
559-571. 

L.S. Da Rios, On the motion of an unbounded fluid with a vortex filament of 
any shape, Rend. Circ. Mat. Palermo 22 (1906) 117. 

G. S. Deem, N. J. Zabusky, Vortex waves: Stationary "V states," interac- 
tions, recurrence, and breaking, Phys. Rev. Lett. 40 (1978) 859-862. 

F. de la Hoz, Self-similar solutions for the 1-D Schrddinger map on the 
Hyperbolic plane, Math. Z. 257 (2007) 61-80. 



33 



[Ho2] F. de la Hoz, Numerical study of a flow of regular planar curves that develop 
singularities at finite time, SIAM J. Appl. Math. 70 (2009), 279-301. 

[HG-CV] F. de la Hoz, C. Garca-Cervera, L. Vega, A numerical study of the self-similar 
solutions of the Schrddinger Map, preprint 2008. 

[Fe] C. Fefferman, Inequalities for strongly singular convolution operators, Acta 

Math. 124 (1970) 936. 

[Fo] D. Foschi, Inhomogeneous Strichartz estimates, J. Hyperbolic Differ. Equ. 
2,1 (2005) 1-24. 

[FM] Y. Fukumoto, T. Miyazaki, Three dimensional distorsions at a vortex fila- 
ment with axial velocity, J. Fluid Mech. 222 (1991) 396-416. 

[GV] J. Ginibre, G. Velo, Scattering theory in the energy space for a class of non- 
linear Schrddinger equations, J. Math. Pures Appl. 9 (1985) 807-836. 

[GP1] R.E. Goldstein, D.M. Petrich, The Korteweg-de Vries hierarchy as dynamics 
of closed curves in the plane, Phys. Rev. Letters. 67 (1991) 3203-3206. 

[GP2] R.E. Goldstein, D.M. Petrich, Soliton's, Euler's Equations, and Vortex Patch 
Dynamics, Phys. Rev. Letters. 69 (1992) 555-558. 

[Grl] A. Griinrock, An improved local well-posedness result for the modified KdV 
equation, Int. Math. Res. Not. 61 (2004) 3287-3308. 

[Gr2] A. Griinrock, Bi- and trilinear Schr dinger estimates in one space dimension 
with applications to cubic NLS and DNLS, Int. Math. Res. Not. 41 (2005), 
2525-2558. 

[GNT] S. Gustafson, K. Nakanishi, T.-P. Tsai, Global dispersive solutions for the 
Gross-Pitaevskii equation in two and three dimensions, Ann. Henri Poincare 
8 (2007), 1303-1331. 

[GRV] S. Gutierrez, J. Rivas, L. Vega, Formation of singularities and self-similar 
vortex motion under the localized induction approximation, Comm. Part. 
Diff. Eqns. 28 (2003) 927-968. 

[GuV] S. Gutierrez, L. Vega, Self-similar solutions of the localized induction approx- 
imation: singularity formation, Nonlinearity 17 (2004) 2091-2136. 

[Har] J.S. Hardin, The velocity field induced by a helical vortex filament, Phys. 
Fluids 25 (1982) 1949-1952. 

[Has] H. Hasimoto, A soliton in a vortex filament, J. Fluid Mech. 51 (1972) 477- 
485. 

[HN] N. Hayashi, P. Naumkin, Domain and range of the modified wave operator 
for Schrddinger equations with critical nonlinearity, Comm. Math. Phys. 
267 (2006) 477-492. 

[HB] E.J. Hopfinger, F.K. Browand, Vortex solitary waves in a rotating, turbulent 
flow, Nature 295, (1981) 393-395. 

[KT] M. Keel, T. Tao, Endpoint Strichartz estimates, Amer. J. Math. 120 (1998) 
955-980. 

[KPV1] C. Kenig, G. Ponce, L. Vega, On the interaction of nearly parallel vortex 
filaments, Comm. Math. Phys. 243 (2003) 471-483. 



34 



[KPV2] C. Kenig, G. Ponce, L. Vega, On the ill-posedness of some canonical non- 
linear dispersive equations, Duke Math. J. 106 (2001) 716-633. 

[Ki] N. Kita, Nonlinear Schrddinger equation with triple 5-functions as initial 

data, Sapporo guest house symposium on mathematics 20 Nonlinear Wave 
Equations (2005). 

[KMD] R. Klein, A. Majda, K. Damodaran, Simplified equations for the interaction 
of nearly parallel vortex filaments, J. Fluid Mech. 288 (1995) 201-248. 

[Ko] N. Koiso, Vortex filament equation and semilinear Schrddinger equation, 
Nonlinear Waves, Hokkaido University Technical Report Series in Mathe- 
matics 43 (1996) 221-226. 

[La] G.L. Lamb, Jr, Solitons on moving space curves, J. Math Phys. 18 (1977) 
1654-1661. 

[LD] M. Lakshmanan, M. Daniel, On the evolution of higher dimensional Heisen- 
berg continuum spin systems, Physics A 107 (1981) 533-552. 

[LF] H. Levy, A. G. Fordsyke, The steady motion and stability of a helical vortex, 
Proc. R. Soc. London. A 120 (1928) 670-690. 

[MHR] T. Maxworthy, E.J. Hopfinger, L.G. Redekopp Wave motion on vortex cores, 
J. Fluid Mech. 151 (1985) 141-165. 

[MB] A. Majda, A Bertozzi, Vorticity and incompressible flow, Cambridge texts 
in Appl. Math., Cambridge University Press (2002). 

[Oz] T. Ozawa, Long range scattering for nonlinear Schrddinger equations in one 
space dimension, Comm. Math. Phys. 139 (1991) 479-493. 

[PeV] G. Perelman, L. Vega, Self similar solutions for modified KdV equation. Some 
applications, J. Differ. Equations 235 (2007) 56-73. 

[PU] F. Planchon, On the Cauchy problem for Besov spaces for a nonlinear 
Schrddinger equation, Comm. Contemp. Math. 2 (2000) 243-254. 

[P12] P. Fabrice, Dispersive estimates and the 2D cubic NLS equation, J. Anal. 
Math. 86 (2002) 319-334. 

[P1V] F. Planchon, L. Vega, Bilinear virial identities and applications, Annal. 
Ecole Norm. Sup., to appear. 

[Sa] P.G. Saffman, Vortex Dynamics, Cambridge University Press, Cambridge 
Monographs on Mechanics and Applied Mathematics, (1992). 

[St] E.M. Stein, Harmonic Analysis, Real Variable Methods, Orthogonality, and 

Oscillatory Integrals, Princeton University Press, Princeton Mathematical 
Series 43 (1993). 

[Ta] T. Tao, A sharp bilinear restriction estimate for paraboloids, Geom. Funct. 
Anal. 13 (2003) 1359-1384. 

[VV] A. Vargas, L. Vega, Global well-posedness for Id non-linear Schrodinger equa- 
tion for data with an infinite L 2 norm, J. Math. Pures Appl. 80 (2001) 
1029-1044. 

[Vi] M.C. Vilela, Inhomogeneous Strichartz estimates for the Schrdinger equation, 

Trans. Amer. Math. Soc. 359 (2007) 2123-2136. 



35 



[We] F.B. Weissler, Asymptotic analysis of an ordinary differential equation and 
non- uniqueness for a semilinear partial differential equation, Arch. Ration. 
Mech. Anal. 91 (1986) 231-245. 

[WD] C. Wexler, A. Dorsey, Contour dynamics, waves, and solitons in the quantum 
Hall effect, Phys. Rev B 60 (1999) 10971-10983. 

[ZHR] N. J. Zabusky, M. H. Hughes, K. V. Roberts, Contour dynamics for the Euler 
equations in two dimensions, J. Comp. Phys. 30 (1979) 96-106. 



36 



